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Abstract 

We investigate azimuthal correlations in deep inelastic diffractive scat- 
tering, e+p— > e + p + X . The dependence of the ep cross section on the 
angle between the lepton plane and some direction in the hadronic final 
state can be written in a simple form; its measurement can be used to 
constrain the cross section for longitudinally polarised photons. Using the 
model of nonperturbative two-gluon exchange of Landshoff and Nachtmann 
we calculate the distribution of the azimuthal jet angle in diffractive dijet 
production and find that useful bounds on the longitudinal cross section 
for such events might be obtained from its measurement. We then discuss 
the predictions of this model for the dependence of the ep cross section on 
the azimuthal angle of the proton remnant p, which contains information 
about the helicity content of the pomeron. 
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1 Introduction 



Our knowledge of diffractive physics in deep inelastic electron-proton scattering 
is greatly increasing as diffractive events are being studied in more and more 
detail at HERA . The phenomenology of these events has many aspects, and 
several theoretical models have been proposed to describe them [|, |], H, ||, ||, . 
Despite various successes of these models we are yet far from a clear theoretical 
picture of what pomeron physics is in terms of QCD. Detailed studies of the 
characteristics of the final state might help to further our understanding of the 
underlying mechanisms and to distinguish between various models. 

The measurement of two different kinds of azimuthal angles has recently been 
proposed: the azimuthal angle of the scattered proton [|, |1^ or, equivalently, of 
the diffractive system as a whole, and the azimuthal angle of the jets in events 
with only two jets of large transverse momentum in the diffractive final state 



1^ . The present paper will be concerned with both issues and has two purposes: 
to discuss some general aspects of azimuthal distributions in diffraction and to 
present in detail predictions for such distributions in the model of nonperturbative 



two-gluon exchange of Landshoff and Nachtmann . 

The structure of this paper is as follows. In sec. ^ we generalise the formalism 
of 1^ for azimuthal distributions in diffraction and show which constraints on the 
7*p cross section for longitudinal photons their measurement can provide. A cor- 
responding framework has long been used in various processes in non-diffractive 
DIS As an application we consider in sec. ^ the azimuthal angle of the 

jets in diffractive dijet production. Some features of its distribution are quite 
characteristic for two-gluon exchange and might offer a way to test the two-gluon 



approximation in this type of events as was already pointed out in ||Tl|. We calcu- 
late the angular dependence in the Landshoff-Nachtmann model and show which 
bounds on the longitudinal cross sections could be obtained from its measure- 
ment. We also show how this method can be generalised to final states that do 
not necessarily have two-jet topology. In sec. |^ we generalise the calculation to 
nonzero t and obtain the corrections this gives for the 7*p cross sections and for 
the distribution of the azimuthal jet angle. Using this calculation we investigate 
a genuine finite-t effect in sec. § the correlation between the azimuthal angles 
of the scattered lepton and proton. In it was shown that this observable con- 
tains information about the helicity structure of the pomeron and argued that it 
might provide a sensitive test of various theoretical ideas about the underlying 
dynamics. We conclude with a summary in sec. |^. 

2 Azimuthal angle dependence in diffraction 

We begin by extending the formalism of to a large class of azimuthal angles 
in diffractive electron-proton or positron-proton collisions. 
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eik)+p{p)^eik')+pip)+X{px) , (2.1) 

where the proton remnant p can be a proton or a diffractive excitation of a 
proton and where four-momenta are indicated in parentheses. We will use the 
conventional kinematic quantities Q^,W^,x,y, s,t for deep inelastic scattering, 
Mx for the invariant mass of the diffractive system X, and the variables /3 = 
gV(g2 + M2 - t) and e = (Q2 + Mi - t)/iW^ + - ml). 

Working in the 7*p rest frame one can write the azimuthal dependence of the 
ep cross section in a simple way by making use of the factorisation of (p.l| ) into 
7* emission by the electron or positron and a diffractive photon-proton collision 

1*{q)+p{p)^X{px)+p{p) . (2.2) 

To achieve this it is essential that the selection of diffractive events, e.g. the 
definition of a rapidity gap between X and p, is unaffected by a common rotation 
about the 7*p axis of the momenta in the hadronic final state Xp, with the lepton 
momenta k and k' being kept fixed. This is guaranteed if the selection criteria 
only involve quantities of the 7*p reaction, i.e. if they do not refer to the lepton 
momenta k and k'. Examples for such criteria have been given in ||^. 

We define an azimuthal angle with respect to a direction in the hadronic final 
state Xp. To this end we introduce a four- vector h which depends only on particle 
momenta in the 7*p reaction ( p. 21) , i.e. on p, q and the momenta of the final state 
hadrons. Using a right-handed Cartesian coordinate system with the z axis in 
the direction of the photon momentum q and some fixed x and y axes we define 
as the azimuthal angle between the lepton momentum k and the vector h, i.e. 
as the azimuthal angle of k minus the azimuthal angle of h. 

We also use h to introduce polarisation vectors for the virtual photon: 

1 ( , Q' \ n 



(2.3) 



where 



{p-q){p-h)-ml {q-h) {p-q){q-h) + Q\p-h) 
hT = h f— q — p 2.4 

is the transverse part of h with respect to p and g, and = e^^''" p^qph^ is 
normal to p, q and h. Polarisations for positive or negative photon helicity are 
as usual given by e± = =f(£^i ± i^2) / V^- 

The contractions of e^, eg, with the appropriate matrix element of the 
hadronic electromagnetic current eJ^ give the amplitudes eM-m for subreaction 
(|2.2|) with photon helicity m, 

eM^m = {Xp out I e J^(0) | p) • , m = -,0,+ . (2.5) 
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The corresponding differential cross sections damm are obtained by multiplying 
•^rn-^m with the phase space element of the hadronic final state Xp and with 
a normalisation factor, summing over the states Xp allowed by our selection 
criteria for the diffractive reaction and averaging over the initial proton spin. 



Our normalisation factor corresponds to Hand's convention |]I5[ for the photon 
flux. From Ai^M..^ with m n we define in an analogous manner differential 
interference terms damn between photons with helicities m and n. It is easy to 
see that the matrix damn is hermitian, damn = da^^. 

With the requirement on the selection cuts formulated above the damn are 
invariant under a common rotation of the momenta in the hadronic final state 
Xp about the 7*p axis. This is because our transverse photon polarisations are 
not fixed but vary with the final state as they depend on h. One can show that 
the cross sections damm are the same for different choices of this vector, whereas 
the interference terms are not. In the following sections we will put extra labels 
on the angle ip and the damn to distinguish different choices of h, though for the 
diagonals damm this would not be necessary. 

Integrating over the phase space of the hadronic final state we obtain a matrix 
Cmn- We will also consider 7*p cross sections and interference terms that are 
differential in some kinematical variables of the final state, such as the momentum 
of the proton remnant or internal variables of the system X. We will only use 
variables that can be defined as Lorentz invariant functions of the four-momenta 
in the 7*p reaction. Provided that the selection criteria for our reaction do not 
refer to any particular frame, the differential cross sections and interference terms 
are then Lorentz invariant and as a consequence depend only on W^, and the 
variables in which they are differential. Due to the rotation invariance property 
just mentioned they are independent of the azimuthal angle of h in our fixed 
coordinate system and hence also of ip. An important property following from 
angular momentum conservation is that interference terms which are differential 
in the direction of h vanish when h becomes coUinear with q and p, in which case 
the azimuthal angle ip is undefined 0. 

The ep cross section can now be written as |^ 

da{ep epX) 2f f 1 , 

-(cr++ + o-__) + ecToo 



dxdQ^dip 27V 12 

— £:cos(2(^) Rea+- +esm{2ip) lma+- 



£{! + £) COS ip Re (cr+o — O"-0 



+ y £(1 + e) sin (y9 Im ((t+q + a-o) 



- TL-Jeil - e) cosip Re (cr+o + cr_o) 
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+ rL'Je{l- e) sin Im (a+o - a_o) 



(2.6) 



where we have integrated over a trivial overall angle, namely the azimuthal angle 
of the scattered lepton in the ep frame, is the helicity of the incoming lepton, 
which is approximated to be massless, e = {1 — y)/{l — y + y^/2) is the usual 
ratio of longitudinal and transverse photon flux and 



2T 



X 



X 



{l-y + y^/2 



(2.7) 



where in the expressions of e and 2r we have neglected terms of order x'^nip/Q'^. 
Equation ( |2.6| ) remains valid if its l.h.s. and the amn on its r.h.s. are made differ- 
ential in additional variables as described above. Since the 7*p cross sections and 
interference terms are independent of if the dependence of the ep cross section 



on this angle is explicitly given by the trigonometric functions in (2.6). 

Let us have a closer look at those combinations of the amn that are multiplied 
with the lepton helicity riin ( |2.6| ). To make their role more apparent we introduce 
differential cross sections and interference terms da^i with k, I = 0,1,2 analogous 
to the damn, but with the linear photon polarisations Eq, Ei, E2 of ( |2.3| ) instead 
of E-, Eq, e+ in the helicity basis. We have the relations 



-{da++ + da—) 

Re (icr+_ 

Im da^ 

Re (licr+o — da-o 
Im (c/o"+o + da-o 



1 



{da J, 



da_ 



Re {da+o + da^o 
Im {da^Q — da^o 



-{dan + da22) 

-]j^{d(Tn - da22) 
Re dai2 

— RerfcTio 
\/2 Re da2o 

— Im dai2 

— \/2 Im da2o 
—\/2\m.daiQ . 



(2. 



The terms that depend on the lepton helicity in the ep cross section are seen to 
be the imaginary parts of 7*p interference terms for linearly polarised photons. 
daki is given by A^^AI; multiplied with a phase space element and real factors, 
summed over the appropriate final states and averaged over the initial proton 
spin, where we define the amplitude eM.k for the reaction 7*^ Xp with photon 
polarisation Ek in analogy to (p75|) . The imaginary parts of daku k 7^ I are 
obviously zero if the phase of this amplitude does not depend on the photon 
polarisation. Note that any (convention dependent) phase of | Xp out) and | p) 
drops out in dan- Also there are no phases coming from eI, Ei because for linearly 
polarised photons the polarisation vectors are purely real; circular polarisation 
introduces extra phases in the 7*p interference terms. 
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We emphasise that in order for Im dan to vanish the do not have to be 
real. The absence of final state interactions gives vanishing interference terms if 
one sums over a set of final states that is invariant under time reversal (cf. ||14||), 
but this is a sufficient condition, not a necessary one. Here we are concerned with 
diffraction and the phases of our amplitudes are certainly nonzero. However, for 
pure pomeron exchange they are given by the signature factor and thus indepen- 
dent of the 7* polarisation. It is beyond the scope of this paper to investigate 
whether for example the superposition of pomeron exchange with exchange of 
secondary trajectories or with multiple exchanges could lead to polarisation de- 
pendent phases that might be tested with longitudinally polarised lepton beams. 

Going back to the for photons with definite helicity, we now make use 
of the parity invariance of strong interactions. It relates amn for different m, n, 
provided that the selection criteria are parity invariant and that h is a vector, 
not a pseudovector. By an argument as in one can show that under these 
conditions 

amniW, Q') = (-1)™+" Q') ■ (2.9) 

Using this and the hermiticity of (Tmn one obtains the relations 

0-++ = (J — , o-+_ = o-_+ = o-|_, cr+o = -a-o , (2.10) 

so that the expression (|2.6| ) of the cross section is simplified: 

da{ep epX) 2T 
dx dQ"^ dip 271 

— 2'^e{l + e) cosv? Rea+o + 2ri^£:(l — e) siny) Imcr+o| . (2.11) 
For linear photon polarisations the relations corresponding to (|2.9|) read 

O"20 = Cr02 = 0-21 = 0-12 = , (2.12) 



(T++ + ecToo - e cos(2v9) 



i.e. transverse photons with polarisation perpendicular to hx of (2]^) do not 



interfere when the final state momenta are integrated over. Expressions analogous 
to ( ^.9| ) to ( f^.l2D are also valid for differential 7*p cross sections and interference 
terms, provided that they depend only on parity even variables, i.e. that one 
sums the damn over a parity invariant set of final states. 



2.1 Bounds on the cross section for longitudinal photons 

We now show how the measurement of the yj-dependence in the ep cross section 
( |2.11|) can be used to constrain the 7*p cross section for longitudinal photons as 
was pointed out in |^. damn is a positive semidefinite matrix, which with the 
simplifications from hermiticity and (|2.9|) from parity invariance implies 

(T++ + (T+_>0, (Too(o-++ - (T+_) > 2|(T+or • (2.13) 
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From the measurement of the (/^-dependence in one can extract the weighted 



sum (Te = cr+^+ecToo of transverse and longitudinal 7*p cross sections as well as the 
interference terms cr+_, Rea+o and Imcr+o- For Imcr+o one needs longitudinally 
polarised electron or positron beams. With unpolarised beams one can use the 
weaker constraints obtained by replacing |(T+o| with Rea+o in ( p.l3[ ) and in the 
following. 



Substituting — eo-Qo for cr++ in (|2.13|) the first condition gives 



^00 < ^ ^ , (2.14) 
whereas the second becomes a quadratic inequality in aoo which leads to 




(2.15) 

Note that the constraints in ( p.l5| ) are the better the larger |o"+o| is compared with 



(Tg — cr+_. This is what one would intuitively expect: a large interference term 
between longitudinal and transverse polarisations implies that neither cxoo nor 
can be very small, i.e. (Too can be neither very small nor very large compared 
with cTg. For |cr+o| <^ cXe — the bounds ( |2.15| ) are less stringent; they then 
differ only by terms of relative order |cr+op/(cre — c'"_|__)^ from the weaker (but 
simpler) conditions 

< o-QO < • (2.16) 



By taking the derivative of ( p.l4[ ) to ( |2.16| ) with respect to e one can see that 



all bounds are decreasing with e if the cross sections (J++, aoo and interference 
terms cr+_, a+o are kept fixed, so that the lower bound is better for smaller e, 
i.e. larger whereas the opposite holds for the upper bounds. Notice however 
that at fixed s a change in ?/ = {W'^ + — vn?^ /{^ ~ f^l) means a change in 

+ and will also change the amn- If their dependence on W"^ + is only 
through a common global factor then this factor drops out in the ratios between 
the bounds on cioo and cioo itself. 

Bounds of the form ( p.l4| ) to (|2.16| ) can also be derived for differential 7*p 



cross sections and interference terms if they satisfy parity constraints analogous 
to ( p.9|) , i.e. if they depend only on parity invariant variables. Thus one can 
obtain bounds on aoo by evaluating inequalities analogous to ( |2.14| ) and ( |2.15| ) 
for differential cross sections and then integrating them. 

The usefulness of the bounds derived here depends of course on how large 
the interference terms are. They will in general be better in some parts of phase 
space than in others, a point we will illustrate in section R^. 
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An important point is that this method allows to constrain the longitudinal 
cross section for fixed e, i.e. fixed y. As is well known, a measurement of the 
longitudinal cross section requires a variation of y, which means that one must 
either measure the ep cross section at different cm. energies ^/s, or, if s is kept 
fixed, have information on how the transverse and longitudinal 7*p cross sections 
depend on + Q^. 

We finally remark that up to now we have not used the requirement of a 
fast outgoing proton or a rapidity gap between X and the proton remnant p in 
reaction (|2.1| ). The analysis developed here, and in particular the possibility to 
constrain the longitudinal cross section, is directly applicable to ordinary deep 
inelastic scattering. In fact, there has been much work on azimuthal correla- 
tions in exclusive or semi-inclusive hadron production and in semi-inclusive jet 



production [13|, with the vector h defining the azimuthal angle (f taken as the 
momentum of the hadron or jet, resp. The (/9-dependence of the ep cross section 
is always given by (|2l6|) , ( |2.11|) with 7*p cross sections and interference terms amn 



appropriately defined for the the process and angle under consideration. 



3 Azimuthal dependence of dijet production 
3.1 Kinematics 

Our first example of an azimuthal angle in diffraction concerns events where 
the proton is scattered elastically and the diffractive final state X consists of a 
quark-antiquark pair at parton level, 

e{k) + p{p) ^ e{k') + p{p) + q{P,) + q{P^) , (3.1) 

which hadronises into two jets. We allow for a finite mass rUg of the quark and 
antiquark. For the vector defining a direction in the dijet system we choose 

h = P = ^iPg-P,) (3.2) 

and work in a reference frame where the incoming p is coUinear with the 7*, 
the photon momentum defining the z axis, and where the total momentum of 
the gg-pair along this axis is zero. It is related to the 7*p cm. by a boost in z 
direction, so that azimuthal angles are the same in both frames. We introduce 
the azimuthal angle (pgg between the electron momentum k and P as in sec. ^ 
and the azimuthal angle x between P and A, where A = p — p is the momentum 
transfer from the proton. We will integrate over x in the present section. Another 
useful variable is the longitudinal component Pl of P, and thus of Pg, along the 
photon direction. Its range is from — p™"^ to P™"^ with pj^"-^ = ^M|/4 — m^, 
and thus independent of x- This is in contrast to the length of the transverse 
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part Pj- of P since 

p2 _ 

1 +t/(M|. -t) -cos^x 



Pt = ^^4/770 7^ V • (3-3) 



Only for t = do the transverse momenta of q and q balance; then is just the 
squared transverse jet momentum. Pl is parity invariant, which will allow us to 
use the simplified expression ( |2.11|) instead of ( |2.6|) to obtain the -dependence 
of the ep cross section. 

Experimentally it is difficult to establish which of the two jets originated in 
the quark q and which in the antiquark q. It is therefore useful to sum over final 
states where the momenta of q and q are interchanged, in other words over P and 
—P. One can show that after this symmetrisation the transverse and longitudinal 
7*p cross sections and the transverse-transverse interference terms are even in P^, 
whereas the transverse-longitudinal interference terms are odd in P^ and vanish 
at Pl = 0. 

To define an azimuthal angle after summing over P and —P one can distin- 
guish the two jets kinematically, e.g. according to which one points in the forward 
direction with respect to the photon. Let Pp be the four-momentum of the for- 
ward and Pb that of the backward jet, and choose for the direction h instead of 

m 

h = PpB = \{PF-PB) (3.4) 

with the corresponding relative azimuthal angle ippB between k and Pfb- The 
longitudinal component of P^b is \Pl\- Writing Pps = sgn(Pi) ■ P we see that 
PpB is a polar vector so that we can again use (|2.11|) for the ep cross section. 



In the case Pl = 0, i.e. when the jet momenta are perpendicular to the 7*p 
axis, eq. ( |3.4|) leaves the sign of h undefined, as the attribution of Pp and Pb 
to the jets is ambiguous. This leads to an ambiguity between azimuthal angles 
(Pfb and ippB + t!' and hence a to sign ambiguity for cos ippB and simppB but not 
for cos{2{p fb) and sin (2^9 ^^5). As mentioned above the transverse-longitudinal 
interference terms vanish at Pl = when summed over P and —P, so that no 
ambiguity remains in the ep cross section. 

The summation over P and —P is in fact trivial under the assumption that the 
diffractive exchange has definite charge conjugation parity, which is of course the 
case for pomeron or two-gluon exchange. Applying charge conjugation invariance 
of the strong interactions to the photon dissociation part of the 7*p subreaction 
it then follows that the differential cross section for remains the same if we 
exchange q and q : 

da {ep ^ ep + q{Pg) + q{Pg)) = da {ep ^ ep + q{Pg) + q{Pg)) . (3.5) 

To sum over P and —P and change from the variable P to Pps we thus only 
need to multiply the cross section with 2 and replace ipgg with ipps and Pl with 
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3.2 The dijet cross section at t = in the LandshofF- 
Nachtmann model 



In this subsection we give the differential cross section for reaction ( |3.1| ) at t = 
in the model of Landshoff and Nachtmann (LN) [l^. The transverse and 
longitudinal 7*p cross sections have been computed in Q], and the calculation of 
the interference terms goes along the same lines. We therefore only recall the 
essentials of the model and give the final results. 

The LN model was developed to give a simple QCD based description of the 
soft pomeron. It approximates pomeron exchange by the exchange of two gluons 
which are taken as nonperturbative, i.e. they have a nonperturbative propagator 
—g^yDip) instead of —g^u/l"^ in Feynman gauge. The nonperturbative propaga- 
tor D{1'^) is a difficult quantity to compute and indeed there is no consensus in the 
literature about its behaviour at small |T^, but we will not attempt to discuss 
this issue here. We will instead follow the rather model independent approach 
of 



12, 16, 18 



which is based on the observation that often the amplitude of a 
considered process can be approximated in such a way that it depends on D{P) 
only via certain simple integrals, so that it is not necessary to know the detailed 
functional form of -D(/^). In the present paper we will need the two moments 

9/3o' 



dl''[af^D{-l^)f 



df[afD{-f)f-f 



47r ' 
Svr 



(3.6) 
(3.7) 



where /3o ~ 2.0 GeV~ and /io ~ 1.1 GeV have been estimated from data [p!6| , p!8 
From the ratio of ( p.7|) and (3^) /io appears as the scale characteristic for the 
behaviour of D onP. Following |TH] we take af^ ~ 1 for the strong coupling in 



the nonperturbative region which dominates the / -integrations in ( |3.6|) , (|3] 

In this model the reaction 7*p qqp is described by the exchange of two 
gluons between a quark in the proton and the gg-pair into which the virtual 
photon splits. In the high-energy limit its amplitude is purely imaginary and 
thus can be calculated by cutting the corresponding Feynman diagrams in the 
s-channel. In each diagram there is then exactly one off-shell quark, namely one 
of the quarks into which the photon dissociates. Its typical virtuality is 



1-P 



{3.t 



which can be seen as the relevant scale of hardness of the process 
obtain the ep cross section from the master equation ( p. 11 

da{ep^epqq) a^m 1 



19 . We 



dipqq dx dQ^ dPi d(3 dt 

, . da'^^_ 
-£ cos(2v9gg-) 



2ttxQ'' 



(1 



dPi dp dt 



2Je{l + e) cosipqg 




+ e 



dPi dp dt 
(3.9) 
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da'^Q/{dPLdp dt) is purely real in our approximation since the Yp amplitude is 
purely imaginary, so that its phase is independent of the photon polarisations, 
see our discussion in sec. ^ For the differential 7*p cross sections and interference 
terms we find 



"''-'■mn 



dPi dp dt 



2 Q^emCq 



(0) 



2(l-ajp(0)) 



(3 



-)2\2 ran i 



(3.10) 



where Cq is the electric charge of the produced quark in units of the positron 
charge and ajpit) = 1 + e + a't with e ^ 0.085 and a' f» 0.25 GeV"^ the soft 
pomeron trajectory. The reduced cross sections 



nqq 



1 -2 



Mi 



Ft + mi 



(Mi + Q^)^L2(P^,w) 



00 



qq 



Mi 



Mi 



5; 



eg? 

'-'+0 



^l + ml 



(Mi + g2)2Li(P^,u;) 



Mi 



P^ + m2 



-2^2 



Q Pl\Pt\ 



M 



X 



Mi 



(Mi + Q^)^Li(P^,u;)L2(P^,w;) , (3.11) 



whose normalisation has been chosen for later convenience, involve loop integrals 



over the functions 



d\l[a^:^^D{-\'^)fMv,w) , 



1,2 



(3.12) 



1 

2w 



v + l-2w 



^{v + 1 - 2wy + Aw{l - w) 
1 



^{v + 1- 210)"^ + 4w(l - w) 



w 



p2 



p2 



(3.13) 



(3.14) 



of the dimensionless variable^ 

Assuming that due to the squared gluon propagator the dominant values of 
in the loop integrals Lj are small compared with we can Taylor expand 
fi{v , w) a.t V = and approximate 

dfiiv,w) 



fi{v,w) ^ V ■ 



dv 



(3.15) 



v=0 



•^The definitions of v and w here differ from those in 
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so that the remaining integral is given by ( p.7| ). The integral ( p.6| ) does not appear 
because both /i and /2 vanish at f = 0. In this approximation becomes 



00 



'-'+0 




1-w) 



l-2w) 



(3.16) 



As a benchmark we have compared the integrals Lj(Py,w) in the approx- 
imation (p. 15 ) with their exact values for the model gluon propagator used in 



le 



D{-P) oc 



1 + 



P 



(n 



n> A 



(3.17) 



where the proportionality constant can easily be obtained from ( |3.6| ). For n — >■ oo 
this becomes D{~P) oc exp(— /^//Xq). We find that the value of n has little 
influence on the Li{P^,w), and that the approximations ( 3.15 ) are in general 
rather good, except however for some regions of parameter space. In particular 
the approximation of Li becomes bad for w close to 1 and for small P|.. On 
the other hand L2 becomes zero and changes its sign for some P|. if w > 1/2 
because the function /2(f) changes sign ai v = w ■ Ix/^t — 2(2w — 1). The 



parameters P^, w for which L2 vanishes are not always well reproduced by the 
approximation; it can be seen that for given P^ the value w = 1/2 obtained 
from (|3.15| ) is too small, so that with fixed the corresponding value of f3 is 
overestimated 

An improved approximation, also leading to the moment ( p.7| ), is achieved by 
replacing (|3.15|) with 



fi{v,w) f» V 



fiivo,w) 



(3.18) 



where we take /q 



With this approximation the values Pj., w where L2 
vanishes are reproduced much better, and the errors on Li are in the region of a 
few percent even if w = 0.9 and 



Pf, as small as 2 GeV^ 
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3.3 Discussion of the results 



Let us make some remarks on the results ( p.l6| ). The first concerns the 

sign of the 7*p interference terms. The transverse-transverse interference is always 
positive, so that the term with cos{2ipqq) in the ep cross section ( p.9|) is negative. 
In |jlT[ it was pointed out that this is the opposite sign than the one obtained 
for gg-production in photon-gluon fusion. The sign of the longitudinal-transverse 
interference depends on the loop integral L2(Pf^, w) and thus on the value of w. 
( |3.16| ) gives a sign change at w = 1/2, the exact value of w from ( |3.11| ) is larger 
and depends on P|. as mentioned at the end of the previous subsection. This 



characteristic change of sign has also been observed in |TT[]. As a general remark 
we can say that the distribution of the ep cross section in (fgg we obtain is very 



similar to the one in the perturbative two-gluon approach of ||TT[ , apart from the 



overall normalisation which comes out different in the two models, cf. [Q. The 
main characteristics of the normalised azimuthal distribution are determined by 
the two-gluon exchange picture. 

We now turn to the dependence on the transverse jet momentum. From 



( |3.11 ) one sees that compared with the transverse cross section the transverse- 
longitudinal interference is suppressed by a factor \Pt\/Mx, the transverse- 
transverse interference by P'^/Mx and the longitudinal cross section by 
(P|n -|- mq)/Mx. This means that the interferences (and the longitudinal cross 
section) are less important if P"^ is small compared with Mx, i.e. if the jets are 
close to the 7*p axis in the reference frame we are working in. Note that for 
light quarks the suppression of the transverse-longitudinal interference is weaker 
than the one for the longitudinal cross section. This suggests a way to experi- 
mentally look for the zero in the longitudinal 7*p amplitude, which is due to the 
behaviour of the integral L2{Pt, w) and can be viewed as a characteristic feature 
of the two-gluon exchange mechanism in this reaction: The zero might be seen 
through the change of sign of the cos ^9^^ -term in the angular dependence as w 
or (3 is varied, whereas it should be difficult to observe it from a dip in the ipqq- 
integrated spectra, given that the longitudinal 7*p cross section is much smaller 
than the transverse one where the zero occurs. We finally note that for heavy 
quarks (|3.11| ) is valid down to = and that all interference terms (though 
not the longitudinal cross section) vanish in this limit where h is coUinear with p 
and q as required by angular momentum conservation. 

For a numerical study we change variables from ipqq and Pl to cppB and \Pl\ 
as explained at the end of sec. ^]l| and integrate over \Pl\- To ensure that we 
have jets and that the scale of ( ^.81 ) remains large we impose a lower cut PTcut 
on P^, which at t = corresponds to an m^-dependent upper cut PLcut on \Pl\, 
see ( |3.3| ). The (/^iT-^- dependence of the ep cross section for the quark fiavours 
u, d, s, c at the HERA cm. energy of ^/s = 296 GeV is shown in fig. |l], where we 
plot da{ep ep qq) / {dippB dx dQ"^ d(3 dt) as a function of ippB for different values 
of the other kinematical variables. We note that in the examples where charmed 
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jets can be produced kinematically their fraction in the 7*p cross section is not 
neghgible; for cases (a) and (6) of the table in fig. |l| it is about 1/3, and for case 
(e) about 1/5. 

The examples in fig. |l| illustrate how a smaller minimum P^/Mj- leads to a 
flatter dependence on ippB as discussed above, while increasing the overall rate. 
The effect of (3 on the sign of the cos ^T-^-term in the cross section can clearly be 
seen. Also shown in the plots is the difference between the approximations (|3.18| ) 
and ( 3.15|) of the integrals in ( 3.11| ): in general the less exact approximation 



( |3.15|) which leads to ( p.l6|) is rather good, especially if 'Pj'cut large. 

Finally we investigate what bounds on the longitudinal 7*p cross section one 
could obtain from measuring the ip i?B-dependence shown in fig. For convenience 
we introduce the quantities 

T^FB _ ""++ T^FB _ ^ "^00 T^FB _ t^FB . t^FB 

^++ " d^dt ' " 'd^dt ' ^ " ^+++^00 , 

whose factors are chosen such that 

dajep^epqq) _ Uem 1 ^ 2 

difpB dx dQ"^ dp dt vr 2'kxQ'^ ^ ' 

{F/^ + F^f cos(2(^fb)+F+^cos(^,.b} . (3.20) 

Up to a global factor F^^ , F^q and F^^ are therefore the Fourier coefficients 
for the v^^^-dependence of the ep cross section, and Fqq^ is the contribution of 
longitudinal photons to Ff^. Table |l] gives the coefficients F^^, F^q and F^^ 
which correspond to the plots in fig. |T], the longitudinal contributions Fqq^, and 
the lower and upper bounds Fj^'^, on F^^ one can obtain from the ^Pfb- 
dependence using the differential analogues of ( ^.15| ). The upper bound ( ^^.14] ) is 



not useful since the transverse-transverse interference da^^/{dp dt) is positive. 

When the minimum P'^/Mx is rather small, i.e. in cases (b) and (e), the lower 
and upper bounds are rather far apart from each other and in this sense not very 
stringent, due to the suppression of the transverse- longitudinal interference by 
\Pt\/Mx compared with the transverse cross section. From ( p.l6| ) we see that the 
lower bound is then suppressed by P'^/M^. On the other hand the longitudinal 
cross section itself has a suppression factor (P^ + m^)/M|^, and as a result the 
lower bound we obtain is quite close to the actual value of Fqq^. 



3.4 Jet angle for more general final states 

We now generalise the jet angle used so far to diffractive final states X that do not 
necessarily have a two-jet topology. As before we work in a reference frame where 
the incoming photon and proton are collinear and where the total momentum of 
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Table 1: Fourier coefficients corresponding to the angular distributions shown in 
fig. |T] and lower and upper bounds -F™, -F^^ on the longitudinal contribution 
Fqq^ one can obtain from them. For the definition of Ff^, F^ff etc. cf. (|3.19|) . 
y = 0.5, e = 0.8, = 296 GeV, t = 0. 







M2 


p2 

' Tcut 








-^00 


pFB 

low 


pFB 

upp 






GeV^ 


GeV^ 






. nb/GeV^ . 






















(a) 


1/3 


80 


16 


-2.9 


-0.75 


4.6 


0.45 


0.22 


0.72 


(b) 


1/3 


80 


4 


-11 


-7.1 


43 


1.3 


0.50 


28 


(c) 


2/3 


20 


4 


-8.3 


3.5 


15 


2.6 


1.0 


3.3 


id) 


1/3 


20 


4 


-54 


-14 


86 


7.2 


4.1 


14 


(e) 


2/3 


20 


1 


-30 


49 


130 


18 


8.4 


80 



X along this axis is zero. Let r be the thrust axis of X in this frame. It can be 
oriented by requiring that it points in the direction of the photon momentum: 

T = Tsgn(q-T) . (3-21) 

This provides a direction in the hadronic final state, which we can also write as 
a four-vector: 

h = (0,t) . (3.22) 
From (|3.21|) and (|3.22|) it follows that /i is a polar vector. Note that in the case of 



a two-jet final state and in the limit t = it becomes proportional to the vector 



PpB defined in (|3.4| ). Another possibility would be to define h from the thrust 
axis in the rest frame of X by equations analogous to ( p.21| ) and ( p. 221 ), and then 
to boost h to the 7*p frame. If X is a dijet this is then proportional to PpB even 
at finite t. 

The vector h defined in one of these ways, or a vector obtained from another 
suitable shape variable of the system X, can be used for the definition of the 
azimuthal angle (f and of 7*p cross sections and interference terms. It is not 
restricted to events with only two jets in X, and it does not require to have jets 
with a transverse momentum large enough for a jet algorithm to be applicable. 
This could allow for a significant gain in statistics. The measurement of the 
(/^-dependence could in particular be used to constrain the cross section for lon- 
gitudinal photons. The discussion in the previous subsection and the numerical 
example with Pf. > 1 GeV^ in table |I] indicates that one might obtain at least a 
useful lower bound even for low Pf^, provided the ratio P'^/Mx is not very small. 
Too small values of P^/M|^ will presumably also present experimental problems, 
since then the polar angle of h is close to zero and the resolution on its azimuth 
will become poor. 
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When the thrust axis is perpendicular to the 7*p direction the requirement 
( |3.21|) does not fix its orientation, so that the angle is only defined up to an 



ambiguity between ^9 and ip + n. This is just as in the case = for the two-jet 
final state which was discussed at the end of sec. |3.1j . Using a similar argument 
as there one can show that the transverse-longitudinal interference terms vanish 
when the thrust axis is perpendicular to q so that no ambiguity appears in the 
ep cross section. 



4 Dijet production at finite t 

4.1 Coupling of the two gluons to the proton 

We will now investigate diffractive production of a jet pair (|3.1| ) at finite t in the 
LN model. Throughout our calculation we take the high energy limit, dropping 
terms that are suppressed by factors of ^. In this approximation t = A'^ where 
At is the transverse part of A with respect to p and q. A characteristic property 
of the LN model is that the two gluons couple to the same quark in a hadron 
[ P^ . The coupling of the gluons to the proton is then given by the isoscalar vector 
current of the nucleon, and the squared amplitude for the process is proportional 

toi 

G\t) = F^it) - ^i^l(t) • (4.1) 

where Fi{t) and F2{t) are the isoscalar Dirac and Pauli form factors of the nucleon, 
respectively, i.e. the sum of the Dirac (Pauli) form factors of the proton and the 
neutron. At t = one has -Fi(O) = 1 and ^2(0) ~ —0.12, cf. 0, and in the region 
\t\ < 1 GeV^ we are interested in Fi(t) is dominating this expression. 

One finds that at high energy the relevant kinematics in the Feynman dia- 
grams for 7*p — i> qqp are determined by the proton momentum and the kinematic 
variables of the 7* — > qq transition, but not on the momentum of the quark within 
the proton. All dependence of the amplitude on the nucleon structure thus comes 
from the form factor G{t) ; there is no further dependence on tranverse or longi- 
tudinal momentum distributions of quarks, even at finite Ay. 

The polarisation vectors for the gluons coupling to the proton both come out 
proportional to the initial proton momentum p. For the complete amplitude they 
are contracted with the propagators of the gluons, and the result is contracted 
with a tensor corresponding to the gluons coupling to the produced quark and 
antiquark. We wish to remark that one need not take a Feynman-like gauge for 
the gluons, i.e. a propagator —g^ijDik'^) where k is the gluon momentum. In 
fact one has some freedom to choose a gauge in our model without changing 
the structure of the result: there is no contribution to the amplitude from the 
tensor k^k^ in a general covariant gauge, nor from fc^rij, + ra^fci^ with some fixed 
four-vector n, which appears in non-covariant gauges. The reason is that in the 
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approximation of our calculation the exchanged gluons couple directly to quarks, 
not to gluons, and thus to a conserved vector current. Note however that for 
terms in the propagator which involve n^riv (appearing in radiative corrections 
to the bare propagator in non-covariant gauges) one would have to investigate 
in detail whether the extra tensors contribute to the leading energy behaviour of 
the amplitude. 

The approximation of two noninteracting gluons in the LN model is certainly a 
very crude one. To go beyond it one could replace the direct coupling of the gluons 
to a quark in the proton by the cut amplitude for the emission by the proton of 
two gluons, in other words by the cut amplitude for g*p — > g*p. Including the 
gluon propagators the latter might be approximated by the gluon distribution in 
the proton [^, ^ at zero t. This is however not useful when we want to compute 
the effects of finite t, since in the gluon distribution the four-momenta of the two 
gluons are by definition equal, in particular they do not transfer any transverse 
momentum Ay. 

Some features of the LN model are also found in this more general framework 
if one makes the assumption that the squared cm. energy of the g*p — > g*p 
amplitude is small compared to W"^ in the region of phase space which dominates 
the amplitude for 7*p qqp. This is for instance the case in the multiperipheral 
approximation. Then one can show that the polarisation of the gluons is again 
proportional to p, and that the relevant kinematics in the Feynman diagrams are 
as in the LN model calculation. Moreover, both statements remain valid if the 
proton dissociates and one integrates over the particle momenta in the proton 
remnant with p being held fixed, provided that p"^ <^ W"^. 

What is however particular to our model is the dependence of the amplitude 
for the emission of two gluons by the proton and their propagation on t and, 
yet more importantly, on the gluon virtualities, the latter being given by the 
nonperturbative gluon propagators. We shall see that precisely these two points 
will have the main effect on the t-dependence of the cross section for our process. 



4.2 Loop integration 

Having contracted the tensor for the two gluons coupling to the proton with 
the one for their coupling to the 7* — > qq transition we must perform a loop 
integration. We label the loop momentum / in such a way that the first gluon 
emitted from the proton carries momentum — / + A/2 and the second gluon 
/ + A/2. Their respective virtualities come out as (/ - A/2)2 = {It - Arf and 
(Z + A/2)^ = {It + At)^, where It is the transverse part of I with respect to p 
and q. Using the cutting rules we are left with a two-dimensional loop integral 
of the form 



L[f] 



(flT 



[affD U\t - AT/2f] D \-{\t + Ar/2)^1 • / 



(4.2) 
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with some complicated function / depending on 1^, Ay and the other kinematical 
variables Py, rrig, Mx, ■ In particular / contains quark propagators whose 
denominators depend on ly ■ and \t ■ Pt so that unlike in the case t = 
we cannot now perform the integration over the angle of It without specifying a 
model for the gluon propagator D. To obtain a more transparent representation of 
the model dependence and to avoid numerical integrations already at amplitude 
level we expand those quark propagators up to second order in (I^ — At/ 2), 
assuming that both and |At| are sufficiently small. The expansion requires 

1^ < A' , llr-Prl < A\ A^ < A\ (4.3) 

where 

^ = _' (4.4) 

1 fj 



is a generalisation to finite t of the scale ( |3.8| ). Note that the first condition is what 
we used in the approximation ( p.l5| ) in sec. p.2| , which was a Taylor expansion 



around If^/A^ = 0. Our calculation will not give an expression analogous to 
(|3.11|) that does not require the gluon virtuality to be small compared with the 
virtuality of the off-shell quark. 

To obtain tractable expressions we also expand denominators in those terms 
which depend on the angle x between P7- and A^. The cross section is then 
a polynomial in cos x and sin x and can easily be integrated over x- The small 
parameter for these expansions is again |Ar|, more precisely they are valid if 

|At-Pt| < + P^ + A|/4 . (4.5) 



4.3 Integrals over the gluon propagators at finite t 



After the expansions just described we are left with a limited number of simple 
loop integrals. They have the form of L[/] in (O) with / = 1, Z^/^^, where 
z,j = 1,2 The corresponding integrands depend only on It and Ay so that the 
integrals are just functions of t. They will turn out to be crucial quantities in the 
discussion of our results in sections and pl2 , 

The integral with / = is zero because its integrand is odd in It, whereas 
the tensor integral with / = /^/^ is related by rotation invariance to integrals 



over the scalars / = I|, and / = (lyAj') . We therefore have three linearly 
independent integrals to evaluate and choose the combinations 



m 
hit) 



L[l] 
L 



L 



2(1t-At)VA| 



It 



1| cos(2(5) 



(4.6) 



where S is the angle between \t and Ay. At t = we have /o(0) = 9Pq/{4:h) and 
/i(0) = 9/9oAto/(8^) froin (P^)> while the integration over 6 gives /2(0) = 0. 
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The ratio of /i(0) and /o(0) involves the scale /Xq ~ 1.2 GeV^, which is also the 
characteristic scale for the t-dependence of the integrals in ( [4.6|) since it is the 
typical scale for the momentum dependence of D{p). We will therefore have two 
kinds of corrections to the cross sections and interference terms at zero t: 

1. corrections in powers of |t| divided by some kinematical variable of the 7*p 
reaction, such as Q^, Mx-, or of ( [4.4|) . By assumption these kinematical 
variables are all large compared with cf. 



2. corrections which depend on t via t/z^o, where /Xq comes from the ratio of 
the integrals in ( [1.6D or from their variation with t. It is important that /Xq 
is not large compared to typical values of |t| < 1 GeV^. 

We have calculated the 7*p cross sections and interference terms keeping the 
corrections in point |l| up to order t/Q^, whereas no expansion was made in t//iQ, 
having in mind that we can have t//io = 0(1). 

One can make a more detailed statement about the small-t behaviour of the 
integrals ( [4.6|) under the assumption that the function Dip) is sufficiently well 
behaved to be Taylor expanded. In (|4.2|) we then can expand (ly =F A-r)^ in the 
gluon propagators around if.. Terms in this expansion that are odd in It-A-t 
vanish after integration over the angle 5, so that the integrals have a power 
expansion in Af. = |t| : 

(1) 1^1 , J2) 




hit) = l + ci^^^ + ci' 



^0 





hit) = ^|t|cf (1 + 4^^^ + ...) . (4.7) 



On one hand ( |4.71 ) shows that the deviation of these integrals from their values 
at t = is proportional to t and not to ^/^. Moreover one may get a reasonable 
description of their t-dependence over a wider range keeping a few terms of this 
expansion. 

We have evaluated the integrals /o(t), hit), hit) with the model ( p.l7| ) of the 



gluon propagator. For n = 4 we obtain good quadratic fits of the integrals in the 
range |t| = to 1.4 GeV^ with 

= -0.5, = 0.12, cS^^ = -0.38, cf^ = 0.09, 

4°^ = 0.027, = -0.31 (4.8) 

and all other coefficients being zero. For n = oo we have an exponential propa- 
gator D(— 1|) oc exp{— l|//io} and easily find /j(t) = exp{ — |t|/(2|Uo)} ■ Jj(0) for 
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i = 0,1,2. Comparing the integrals for n = 4 and n = oo we find that they 
are almost equal for Jq and that Ji is smaller for n = oo than for n = 4, the 
coefficient for n = oo being —0.5 instead of —0.38. I2 vanishes at all t for the 
exponential propagator, for tt, = 4 it is still very small compared with \t\ ■ Iq. 



4.4 Results 

We now present our results for the 7*p cross sections and interference terms 
of quark-antiquark production at finite t, the ep cross section is obtained from 
eq. (|3.9|). We first give analytical expressions including the corrections in t/ii^, 
but without the corrections of order t/Q"^ which are rather lengthy. The latter 
will be included in the numerical discussion below. 

To zeroth order in t/Q"^ only two linear combinations of the integrals ( |4.6| ) 
appear, namely 

Ki{t) = hit) - \t\ /o(t)/4 , K2{t) = hit) - \t\ hit)/ A . (4.9) 

From (^4.7|) we see that the leading term in the expansion of Ki it) in t//io is 
constant, whereas the leading term for K2 is proportional to t//io- We introduce 
the abbreviations 



Ml/A-t/A- Pi Vl + ml 

(4.10) 



M\ - 1 t^o M\ 
^ Ml/ A -ml- Pi VI 

Ml/A - t/A -Pl^ Vl + ml ' 

whose limits for t ^ are given for easy comparison with our results in sec 
and the variable = (1 — /5) 6 as a generalisation to finite t oiw defined in (|3.14]) . 
The result then reads 

van 

dPLd/3dt ~ 

a^^el ^ G\t) , ^ ^ SZ (4.11) 

3 (3 J Ml - tiMl + Q2 _ t)2 



with 



S 
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4(l-2a)6(l-w;)2 + (l-6)(l-2w)' 



a2 



+ 



K 



a 



[1 - 2a) b il - 2w + 2w^) + 2il - b)w^ + O(VQ^) , 



S'ol = -j^^\{AKlil-2w)' + mW} + Oit/Q'^ 



Sf_ = -{SKlhil-w)^ -AKlhwil-w)] + Oit/Q^ 
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[kI (1 - - 2w) - Kl w{l - 2w)/2] + 0(t/g^) . (4.12) 

For t ^ we recover our previous expressions (|3.1CI|) , (|3.16|) . 

Let us now give some numerical examples, obtained with the parametrisation 
([4.8| ) of the integrals /j(t) which corresponds to the model gluon propagator 
(|3.17|) with n = 4. As in sec. ^]3| we change variables from Pi and ipqq to \Pl\ 



and LfFB and integrate over \Pl\- We impose an upper cutoff on \Pl\ and sum 
over the three light quark flavours m, d, s. In fig. |^ we plot the t-dependence of 
the quantities Ff^, F§P , Fff , F^^ introduced in ( p7T9|) , but taking out the 
squared proton form factor G^(t) and the t-dependent part of the Regge 

power, both of which give a rather strong suppression of the cross section at t 
away from zero. As a result of the different t-behaviour of F^^ , F^^ and F^^ 
the (/JiT-B-dependence of the ep cross section will change with t. 

In order to see to what extent these results depend on the specific form of the 
gluon propagator we plot in fig. ^ the same quantities as in fig. 0, now with the 
simplest ansatz for the integrals we can make: 

= AW = ^f. /.W = ^Wcf. (4.13) 



keeping only the lowest order in t of the expansions ( [4.7|) . The leading coefficient 
€2'^ in I2 is not determined from phenomenology as is the case for Iq and Ji, and 
we take three different values 0, 0.5 and —0.5. We see how the variation of 0^2^ 
modifies the behaviour of the Fourier coefficients at moderate and large values of 
\t\ quite drastically; it can for instance lead to a change of sign in the interference 
terms at fixed t. One would however have to see whether such large variations of 
€2'^ can be obtained with realistic gluon propagators. 

Comparing the plot for C2°^ = and the corresponding one obtained with our 
special ansatz for the gluon propagator for which !^ we see that the effect 
of approximating the t-dependence of the integrals Jj by the leading terms in the 
expansions (|4.7| ) is by no means small. This is not surprising as the leading order 
approximation is only expected to be good for \t\/fil <^ 1. Notice also that the 
first order coefficients Cq^ and Ci^^ in ([4. 8]) are rather large. 



We have compared the results which include corrections up to order t/Q'^ with 
the expressions given in ( [4.12| ) where only the t-dependence through t/fiQ in the 
integrals Ki(t) is kept. In most of parameter space the latter give a very good 
approximation, and even for rather small Q^, or rather low minimum 
the formulae ( ^.12| ) give the correct qualitative features. Apart of course from 
the squared elastic form factor G{tY and the t-dependent pomeron trajectory the 
main effect in the t-dependence of the 7*p cross sections and interference terms 
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thus turns out to be from the integrals Ii{t), i.e. from the fact that at t 7^ the 
two exchanged gluons have different virtuahties. On one hand this means that 
corrections in t/Q^ are less important in the kinematical region we are investigat- 
ing. On the other hand the results depend on the details of the nonperturbative 
gluon propagator encapsulated in the Ii{t), and the phenomenological constraints 
(|3.6| ), ( p.7| ) are not sufficient to predict the t-dependence quantitatively, they only 
provide the right order of magnitude and the characteristic scale /Xg. 



5 The azimuthal angle of the scattered proton 
or proton remnant 

In this section we turn our attention to another azimuthal angle in diffractive 
processes choosing for the vector h 



h 



Px 



(5.1) 



(fx is the azimuthal angle between the lepton and the diffractive system X, i.e. 
ipx + vr is the azimuthal angle between the lepton and the scattered proton or 
proton remnant. This angle was introduced and discussed in 0. Note that the 
7*p cross sections and interference terms a^J introduced there are integrated 
over the internal momenta of the system X but not over p. In the notation used 
in this paper they read 



a, 



1 f3 da, 



X 

mn 



rcf . 



vr ^ d(3 dt 



(5.2) 



this paper 



Px becomes collinear with q and p when \t\ takes its minimum value, which 
is zero in the high energy limit. As we remarked in sec. |^ the 7*p interference 
terms must then vanish for t — because of angular momentum conservation, 
and the crucial question we will be concerned with in the following is how fast 
they do. To quantify this we normalise the interference terms with respect to the 
7*p cross sections and consider the ratios 



R+o — 



da^_/{dpdt) 
da^+/{d(3dt) + da^o/{d(3dt) 

da^Q/{dpdt) 
da^+/{d(3dt) + da§,/{d(3dt) 



(5.3) 



If they behave like a (possibly fractional) power of |t| for t — > then the scale 
that compensates t in these dimensionless quantities determines how large they 
are at finite t. 
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It was argued in and confirmed by an explicit calculation in that in 
the phenomenological pomeron model of Donnachie and Landshoff pS) one has 



fl.- ~ ;P , fi« ~ if . (5-4) 

where Q could be replaced by Mx or some combination of Mx and Q, the im- 
portant point is that the scale dividing t is a kinematical quantity of the 7*p 
subreaction, and therefore rather large compared with t for the typical values of 
t, M|- and in diffractive DIS. 

For models that describe diffraction in terms of soft colour interactions [[7|, 
or in models where the QCD vacuum plays an important role [53], one can expect 



a different behaviour . In these models there is some scale characteristic of soft 
physics which could take the place of Q in the expressions of (|5.4|). This would 
lead to larger interference terms and thus to a more pronounced (/3x-dependence 
of the ep cross section. In the LN model the diffractive mechanism is described 



by soft gluon exchange; we will see in sec. |5.2| where and when its typical scale 
/io replaces Q in ( |5.4| ). We remark that the powers of \t\ which give the small-t 
behaviour of and -R+o niay in general be different from those in ( |5.4| ). 

5.1 Calculation in the LN model 

We now turn to the predictions of the LN model for the dependence of the ep cross 
section on the proton angle. We first have to replace the general diffractive final 
state X with a quark-antiquark pair. This is the lowest order approximation of X 
at parton level and should give a reasonable description, except in the region of 
small /5, or large diffractive mass Mx, where additional gluon emission is known 
to be important. 

The calculation of the 7*p cross sections and interference terms for the process 
( |3.1|) is essentially the same as the one in sec. | with P of eq. ( |3.2|) replaced by px 
in the expressions of the photon polarisation vectors ( p.3|) . We integrate again 
over the relative azimuthal angle between P and px , but now the azimuthal angle 
between k and px is kept fixed, not the one between k and P. 

We have to make an additional restriction on the diffractive final state, be- 
cause we need that of eq. (|3.8|) , ( |4.4|) is large for our approximations described 
in sec. O to be valid. Unless we have a large mass for the produced quarks. 



this means that we must impose a lower cut on their transverse momentum. 

In [0 it was shown that the 7*p cross sections and interferences have a physical 
interpretation in terms of the helicity of the pomeron if one works in the rest frame 
of the diffractive system X, provided that the selection cuts on the hadronic 
final state are invariant when the particle momenta in the system X are rotated 
around the photon momentum in this frame while all other momenta are kept 
fixed. Nonzero interference terms between photons with definite helicities in 
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this frame imply that different amounts of angular momentum along the photon 
direction are transferred from the proton. In pomeron language this means that 
the pomeron can carry different helicities. 

To satisfy the above criterion of rotation invariance we impose a cut on the 
transverse quark momentum in the X rest frame and not in the 7*p system. To 



do this we have to transform the kinematical quantities introduced in sec. ^ 
to the cm. of X. We denote three-momenta with an asterisk there and use a 
right-handed coordinate system with the z axis along the photon momentum q*. 
Both the transverse and the longitudinal momenta of the gg-pair are opposite to 
each other, not only the longitudinal ones as in the frame we used in sec. 
Thus P* = (P* — P|)/2 is equal to the three-momentum of the quark jet. Instead 
of X ci'iid Pl of sec. we use in the X rest system the relative azimuthal angle 
X* between P* and p*, and the longitudinal momentum P£ of P* = P*. Its 
kinematical limits are the same as for P^. For the transverse component P^ of 
P* in this system we have P^^ = M\/A — — P£^, compared with (^31). The 
relation between the longitudinal and transverse components of P in both frames 
is as follows: 

Pl = . ^ i^^ (P2-|P^|cosx*-(l-2/j)v^/Mx) , 



|Pt|cosx = - ; ^ :(|P^|cosx* + P2-(l-2/3)v^/Mx) , 



|Pr|sinx = — |P^|sinx* . (5.5) 

Integrating over x* we obtain the ep cross section differential in y^x, a;, Q^, t, /3 
and P£. Note that P£ and x* are defined in the X rest frame, but in the 7*p 
system. 



5.2 Results 

From our master formula the (y9x-dependence of the ep cross section is given 

by the analogue of ( p.9|) in sec. with the replacements ipx^ Pl — ^ Pl 

and da'^^ da^^. The transverse-longitudinal interferences are real for the 
same reason as in the case of the jet angle. As in sec. |4.3| we have calculated the 
7*p cross sections and interference terms up to order t/Q"^, treating t//ig as of 
order 1. Again we will not give the analytic expressions of the 0{t/Q^) terms, 
but use them in our numerical discussion. We find 

""^ mn 



dPld(3dt 
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with 



S 



X 



n *2 



,*\2 



'^00 



4(1 - 2a*) - w*Y + (1 - h*){l - 2w 
(1 - 2a*) h*{l - 2w* + 2w*'^) + 2(1 - h*)w*'^ 
1 (1 - 2w*)' + ml w*'} + 0{tlQ 



+ 



0(t/Q2 



5: 



5 



X 



8KIK2 

— L^6* 1 

a* ^ 



,*\2 



+ 0{t/Q' 



1 



*2\ 



Ml - t a*2 

-(l-w*)(l-7w* + 8w;* , 
+4a*(l - w*){l - 8w* + lOw*^) 
+2(1 - 4a*)(l - /3)(1 - &*)(1 - 6m;* + 6^*^^ 



+K1K2 [-{1 - w*){l - llw* + 16w*2)/2 
+2a*(l - w*){l - 12w* + 20w*^) 
+ (1 - 4a*)(l - (3){1 - b*){l - lOw* + 12w*^) 

+Kl \w*(2 - 9w* + 8w*2)/2 - Aa*w*{l - 5w* + bw*"^ 
-2(1 - 4a*)(l - /3)(1 - 6*)w*(l - 3w*) ] } 



Here we have used the integrals -ft'i(t), -ft'2(i':) introduced in 
viations 



(5.7) 
and the abbre- 



TD*2 



p i2 _L ^2 



+ m 



P5? + m 



, *2 



;i - /3) 6* 



Mi 



(5.8) 



A remark is in order on the appearance of the square root y— t in the ex- 
pression of the transverse-longitudinal interference terms. One might suspect 
that there is a contradiction with the analyticity properties of scattering ampli- 
tudes, but this is not so. The point is that this interference term is multiplied 
with cos ipx in the ep cross section, and that the expression of cos ipx in terms 
of Mandelstam invariants also involves square roots, so that the appearance of 
is a consequence of the kinematical variables we choose. Put in a different 
way, the interference terms can have a dependence on through the polari- 
sation vectors ( p.3| ), which contain square roots. A corresponding remark can be 



25 



made for the appearance of Pl in the transverse-longitudinal interference term 
corresponding to the jet angle in (|3.11|) , ([4.12|) . 



As discussed in sections |4.3| and 14.41 we have in the hmit of small t 



where c^^ is not known from phenomenology and has to be obtained using a 



specific ansatz for the nonperturbative gluon propagator. With ( |5.7|) we find 
for the small-t behaviour of the interference terms normalised to the 7*p cross 
sections 



~ (cf - 1/4) 4 , i?+o ~ ^ • (5.10) 

Note that the behaviour of in the limit t — is determined by the coefficient 
\ in contrast to the case of the jet angle investigated in sec. [4.4| , where we 
could make a parameter-free prediction for this limit. In particular the sign of the 
transverse-transverse interference depends on the details of the gluon propagator. 
If is close to 1/4 one even finds that the 0{t/Q^) terms are dominant and 
-R+_ ~ Apart from this caveat the result (|5.10| ) is however independent 

of the detailed properties of the gluon propagator. The scale /Xq comes into play 
in this model via the nonperturbative dynamics of the exchanged gluons. Since 
it only appears squared in the calculation, cf. ( [4.7| ), it is clear that it can not be 
the scale dividing in -R+o; there we have the large kinematical variable Q as 
in Q. 

The situation is however more complicated for than discussed so far. 
We will not find ratios -R+_ of order one as suggested by ( ^.10|) in our numerical 



examples. The reason is that (but not S^q) has an additional suppression 
compared with the cross section term S^_^^ by a factor of a*b* = P^/Mx, which 
can be rather small for the /M^ where we cut in the case of light quarks, and 
even goes down to zero for charm production. After integration over P^^ the first 
relation of ( |5.10| ) is more precisely 



R^- - (4°^ - 1/4) 1^ ■ ^ , (5.11) 



for light quarks with a momentum cutoff and 

~ (4°' - 1/4) § . ^ . (5.12) 

for heavy quarks and integration down to = 0. In (|5.12| ) appears since it 
is the typical scale of P^-^ in the integration if the quark mass is large. 

The question arises what one can expect for when there is no lower cut on 
P^ in the case of light quarks. After all the main contribution to da^^/ {dt df3) 
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summed over all flavours and the full phase space is from light quarks at low 
transverse momenta. Because of the approximations of our calculation we cannot 
extrapolate ( |5.7| ) to this region, but we want to give an educated guess. We have 
argued in [|], that with some caveats the LN model can still be applied to 
7*p qqp in the limit + m'^ — >■ 0. In the results ( |3.11| ) from our investigation 
of the jet angle dependence at t = we observe that the suppression of S'^_ with 
respect to 1S4.+ is by a factor P'^/Mx even in this limit. Taking this as a guidance 
for the interference term in our present problem we expect that a suppression by 
P^ /Mx of the differential interference term may persist for very small P^-^ + 
so that Ptcui ( |5.11|) is to be replaced with some average P^? if we integrate 



over the full phase space. Examining the loop integrals LiiP^^w) of ( ^.12[ ) one 
further finds that the typical scale for the P\ -dependence of the cross section is 
/Xq in the case where is small and (3 not too close to 1. This leads us to the 
guess 

for the interference term without a cut on P^^, which is the quantity originally 
discussed in 0. Notice that the Oit/Q"^) terms in now also contribute 
to the leading term of -R+_. The scale /Xq has cancelled and we find the same 
behaviour for both and i?+o as in the Donnachie Landshoff model (|5.4|) . Let 
us however remark that there one has ~ \t\/Q'^ even with a large cutoff on 
P^^, in contrast to ( p. 101) , (|5.11|) , so that the predictions of the two models are 
by no means identical. 

Coming back to what we were able to calculate in the LN model we now give 
some numerical illustrations of our results. In analogy to sec. iA we integrate over 
PI and plot the Fourier coefficients , F^_, F^q with a global factor &{t) 
taken out. They are defined like Ff^, F^^, in (|3.19|) with the superscript 
FB replaced by X and thus appear in da{ep epqq)/{dipx dx dQ"^ dj3 dt) in a 
manner analogous to ( p. 201 ). We either sum over the three light fiavours u,d,s 
for the produced quarks, with a minimum P^ so that our calculation is valid, 
or we consider produced charm quarks for which we can integrate over the full 
kinematical range of P£. 

The t-dependence of the Fourier coefficients for the model propagator ( |3.17| ) 
with n = 4 is shown in fig. ^for different values of the free kinematical parameters. 
We observe that the transverse-transverse interference is usually larger than the 
transverse-longitudinal one. Fig. ^ shows an example of the yjx-dependence of 
the ep cross section at two different values of t. In case (6) the distribution is 
clearly not flat, although the effect is not very large, whereas in case (a) almost 
no v^x-dependence can be seen. This illustrates how the transverse-transverse 
interference is affected by the parameter P^^^^/Mx, which is 1/20 in case (a) 
and 1/5 in case (6). 

To assess the model dependence of our prediction we also evaluated the Fourier 
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coefficients taking the simple ansatz ( [4.13| ) for the integrals over the gluon prop- 
agators, with different values for the coefficient C2^. They are shown in fig. ^ 
for the case of light quarks, the effects for charm are similar. As in sec. the 
results, especially at large \t\, change considerably with C2°^ In particular the 
sign of the transverse-transverse interference term is different for below or 
above 1/4, and for cf^ = 1/4 this term is very small, as discussed above. We 
repeat that one would have to see whether realistic gluon propagators give values 
of as far away from zero as the ones taken in fig. ^. 

We find that the leading order expressions of <S^^, Sqq and in ( |5.7| ) give 
a rather good approximation of what is obtained by including terms of 0{t/Q'^), 
except of course for if c"2^ = 1/4. As in sec. [4.4| this means that the main effect 
comes from terms depending on t/fi^, whereas corrections in t/Q"^ are relatively 
small. Terms in t//iQ are also essential to describe the t-dependence of S^q. 
Although its order of magnitude at small t is given by v^^/Q, a mere square 
root dependence on |t| for the Fourier coefficient F^q in figures ^ and ^ is clearly 
not a good approximation unless t is very small. 

To conclude this section we remark on the possibility to constrain the cross 
section for longitudinal photons from the measurement of the (yJx-dependence 



using the method described in sec. 2.1. With the results we obtain in the LN 



model the bounds on daoo/{dt d(3) would not be stringent at all, and be far away 
from its actual value. This is because we find the interference da+o/{dt df3), whose 
size is crucial to obtain good constraints, to be of order -y/jtl/Q. From ( |2.16| ) we 
can see that the lower bound on daoQ/{dt df3) then vanishes like {tl/Q"^ for small t, 
whereas daoo/{dt d(3) itself does not become small in this limit. This is different 
from the situation we found for the jet angle in sec. |3]^. One would expect better 
bounds if the ratio of the longitudinal-transverse interference and the 7*p cross 
sections at small t were dominated by a hadronic scale instead of Q or Mx- This 
might happen in other models of diffiaction where soft dynamics is important. 



6 Summary 

In this paper we investigated correlations between azimuthal angles in deep in- 
elastic ep diffiaction, using the one-photon approximation. We ffist derived the 
general expression for the dependence of the ep cross section on a suitably de- 
fined azimuthal angle between the lepton plane and a direction in the hadronic 
final state in terms of cross sections and interference terms of the 7*p collision 
for different photon helicities. This was a direct generalisation of the work in 
We showed that those terms in the cross section that depend on the helicity of 
the lepton beam are sensitive to a polarisation dependence of the phases in the 
7*p amplitudes for linearly polarised photons. From the angular dependence of 
the ep cross section one can obtain bounds on the differential or integrated cross 
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section for longitudinal photons, without having to vary y as it is needed for its 
direct measurement. How stringent these bounds are depends on the size of the 
interference term between longitudinal and transverse polarisations and thus on 
the choice of azimuthal angle and on the region of phase space considered. 

We have investigated the dependence on the azimuthal jet angle predicted 
by the LN model for the parton level reaction ep — > ep + qq at large transverse 
momentum of the gg-pair, which at hadron level describes a pair of jets that 
carries the entire four-momentum of the diffractive final state. The size of the 
interference terms is found to be controlled by the quantity P|./M|^. The sign 
of the transverse-longitudinal interference depends on (3. Since this interference 
is less strongly suppressed than the longitudinal cross section, it may offer an 
opportunity to observe the zero of the longitudinal amplitude at certain values 
of and (3 which is characteristic of the two-gluon exchange mechanism. The 
bounds on the longitudinal cross section obtained from the azimuthal dependence 
might be quite useful, at small Pf,/M|^ at least the lower bound comes out quite 
close to its actual value which is also small in this kinematical region. We suggest 
that the use of an azimuthal angle defined from an event shape variable like the 
thrust axis in the diffractive final state would allow to extend this method to a 
wider class of final states, in particular it would allow to go to smaller values of 
Pf. than those needed for jet algorithms and thus to increase the total rate in the 
analysis. 

The cross section for ep ^ ep + qq was then calculated at finite t with the 
approximations A|, < and \\ < A^, for the definition of A^ cf. (^^. 
Its region of validity is therefore the production of jets or heavy flavours where 
A^ is sufficiently large. The result then involves three t-dependent integrals ( [4.(j| ) 
with two gluon propagators at different virtualities. The relevant scale for the 
t-behaviour of these integrals is /Iq ~ 1.2 GeV^. The limit t — > for two of them 
is known from phenomenology, for the rest one has to resort to specific model 
propagators. 

Applying this calculation to the dependence on the azimuthal jet angle we 
find that apart from the dominating effect of the proton form factor G{t) and 
the pomeron trajectory ajp(t) the t-dependence of the 7*p cross sections and 
interference terms is controlled by corrections in t//iQ coming from the gluon 
propagators, corrections in t divided by a large kinematical scale of the transition 
7* — i> qq are much smaller. As a consequence the quantitative features of the 
results depend on the choice of gluon propagator. Using the model propagator 
(|3.17| ) we typically find that the sum of the transverse and longitudinal cross 
sections decreases by a factor around 2 between |t| = and 1.4 GeV^ when the 
strong suppression from G{t) and the pomeron trajectory is taken out. The 
absolute size of the interference terms tends to decrease with \t\ and one can even 
have a change of their signs. 

Another important azimuthal angle is that of the scattered proton or proton 
remnant. It was shown in that its measurement can give information on the 
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helicity structure of the pomeron. We have investigated its distribution in the 
LN model, but due to our approximations had to restrict ourselves to a qq final 
state with large transverse quark momentum P^^ in the qq rest fi:ame, or with 
a large quark mass. Like for the jet angle we find that the t-dependence of the 
cross sections and interference terms is controlled by the scale /Iq, and that the 
results depend rather strongly on the integrals over the gluon propagators. The 
order of magnitude of the transverse- longitudinal interference is given by y/^/Q. 
The ratio between the interference of the two transverse polarisations and the 
transverse cross section goes like t/fi^ which can be large, but it is suppressed 
by an additional factor P^/M]^ so that this interference is small at low P^^. 
Unfortunately we cannot take the limit P^^ — > for light quarks in our calculation 
but our guess is that the P^? -integrated transverse-transverse interference will 
be suppressed by t/Q"^ compared with the transverse cross section, which would 
lead to a rather fiat angular dependence in the ep spectrum. According to the 
discussion in |^ the helicity of the LN pomeron is then dominated by one value 
in the inclusive diffractive process, whereas several helicities are important when 
there is a high transverse momentum or mass scale in the diffractive final state. 
Notice that a system of two gluons we use to model the pomeron can in principle 
transfer any integer value of angular momentum through its orbital motion. 

Our finding that finite-t effects are rather sensitive to the nonperturbative 
gluon dynamics in the LN model suggests that they may come out quite different 
in other models of diffraction and could thus be a useful probe of the mechanisms 
at work in diffractive physics. 
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Figure 1: Dependence on ippB of da{ep epqq) / {dippB dx dQ'^ d(3 dt), summed 
for M, d, s and c quarks. Values of the free kinematical variables are ^/s = 
296 GeV, y = 0.5, e = 0.8, t = and those given in the table. A cut ^ < 0.05 
has been imposed. Full lines correspond to the improved approximation (|3.18| ) 
of ( |3.11| ), dotted ones to the simplified results ( p.l6| ) obtained with the approxi- 
mation ( p. 151) . 
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Figure 2: Fourier coefficients F^^ , , F^^ in the ep cross section and the 
contribution Fqq^ of longitudinal photons to F^^ . For their definition cf. ( 3.19| ). 
They are summed for u, ci, s quarks and a global factor * is taken out in 

the plot. The results here are obtained with the model gluon propagator (|3.171 ) 
for n = A. Kinematical variables are (a): ^ = 296 GeV, y = 0.5, = 40 GeV^, 
/3 = 1/3, \Pl\ < 4 GeV. (6): y/s and y as before and = 40 GeV^ p = 2/3, 
\Pl\ < 1 GeV. 
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Figure 3: As fig. ^ (a) but with the ansatz ( |4.13| ) for the integrals over the gluon 
propagators with different values of Remember that a factor &{t) is 
taken out in the plot; the differential cross section does not rise with |t|. 
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in the ep cross section, defined in 
analogy to ( ^.191 ), with a global factor &{t)^~'^"'^ taken out. Note that F^ is 
scaled down by a factor 10. (a) and {b) are summed for u, d, s quarks with a 
lower cutoff P^^ > 4 GeV^, (c) and (d) are for charm quarks without a cut on 
P^. Kinematical variables are ^/s = 296 GeV, y = 0.5 in all cases and (a): 
g2 = 40 GeV^ /3 = 1/3. (6): = 40 GeV^ /3 = 2/3. (c): = 25 GeV^ 
/3 = 1/3. (d): = 6.25 GeV^, P = 1/3. The curves are obtained with the 
model gluon propagator (|3.17| ) with n = A. 
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(6): The same for the Fourier coefficients from fig. ^ (6) 
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Figure 6: As fig. § (6), but with the ansatz ( 4.13| ) for the integrals over the gluon 
propagators with different values of cfK 
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